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A qualitative analysis of the chiral phase transition in QCD with two massless quarks and non-zero baryon 
density is performed. It is assumed that at zero baryonic density, p = 0, the temperature phase transition 
is of the second order. Due to a specific power dependence of baryon masses on the chiral condensate the 
phase transition becomes of the first order at the temperature T = T p h(p) for p > 0. At temperatures 
Tcont(p) > T > Tph(p) there is a mixed phase consisting of the quark phase (stable) and the hadron phase 
(unstable). At the temperature T = T con t(p) the system experiences a continuous transition to the pure 
chirally symmetric phase. 
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It is well known, that the chiral symmetry is valid in perturbative quantum chromodynamics (QCD) with 
massless quarks. It is expected also, that the chiral symmetry takes place in full-perturbative and nonperturba- 
.tive QCD at high temperatures, (T > 200 MeV), if heavy quarks (c, 6, t) are ignored. The chiral symmetry is 
lO strongly violated, however, in hadronic matter, i.e. in QCD at T = and low density. What is the order of phase 
-transition between two phases of QCD with broken and restored chiral symmetry at variation of temperature 
and density is not completely clear now. There are different opinions about this subject (for a detailed review 
see Ref. [fl |2] and references therein) . 

In this paper we discuss the phase transitions in QCD with two massless quarks, u and d. Many lattice 
i calculations 0103121 indicate, that at zero chemical potential the phase transition is of the second order. It will 
be shown below, that the account of baryon density drastically changes the situation and the transition becomes 
r*| pf the first order, and, at high density, the matter is always in the chirally symmetric phase. 

Let us first consider a case of the zero baryonic density and suppose that the phase transition from chirality 
violating phase to the chirality conserving one is of the second order. The second order phase transition is, gen- 
erally, characterized by the order parameter 77. The order parameter is a thermal average of some operator which 
-may be chosen in various ways. The physical results are independent on the choice of the order parameter. In 
QCD the quark condensate, 77 = |(0|m/|0)| = |(0|dd|0)| ^ 0, may be taken as such parameter. In the confinement 
phase the quark condensate is non-zero while in the deconfinement phase it is vanishing. 

The quark condensate has the desired properties: as it was demonstrated in the chiral effective theory pHIEj, 
77 decreases with the temperature increasing and an extrapolation of the curve 77 (T) to the higher temperatures 
indicates, that 77 vanishes at T — « 180 MeV. Here the superscript "0" indicates that the critical tempera- 
ture is taken at zero baryon density. The same conclusion follows from the lattice calculations [HI E3 Ej , where it 
was also found that the chiral condensate 77 decreases with increase of the chemical potential [TU1 ITT] . 

Apply the general theory of the second order phase transitions 12 and consider the thermodynamical po- 
tential $(77) at the temperature T near Ti°\ Since 77 is small in this domain, $(77) may be expanded in 77: 

*(V) = $0 + \a 77 2 + \b t/ 4 , B > . (1) 

For a moment we neglect possible derivative terms in the potential. 

The terms, proportional to 77 and 77 s vanish for general reasons I n QCD with massless quarks the absence 
of 77 and 77 s terms can be proved for any perturbative Feynman diagrams. At small t = T-T c (0) the function A(t) 
is linear in t: A(t) = at, a > 0. If t < the thermodynamical potential $(77) is minimal at 77 7^ 0, while at t > 
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the chiral condensate vanishes 77 = 0. At small t the t- dependence of the coefficient B(t) is inessential and may 
be neglected. The minimum, 77, of the thermodynamical potential can be found from the condition, d^/dr] = 0: 



{ 0, * > 0. 

It corresponds to the second order phase transition since the potential is quartic in 77 and - if the derivative terms 
are included in the expansion - the correlation length becomes infinite at T = . 

Turn now to the case of the finite, but small baryon density p (by p we mean here the sum of bar yon and 
anti-baryon densities). For a moment, consider only one type of baryons, i.e. the nucleon. The temperature of 
the phase transition, T p h, is, in general, dependent on the baryon density, T p h = T p h(p), with T p h(p = 0) = 
At T < T p h(p) the term, proportional to Ep, where E = ^ p 2 + m 2 is the baryon energy, must be added to the 
thermodynamical potential QJ. As was shown in ^3 the nucleon mass m (as well as the masses of other 
baryons) arises due to the spontaneous violation of the chiral symmetry and is approximately proportional to 
the cubic root of the quark condensate: m = cr] 1 / 3 , with c = (87T 2 ) 1 / 3 for a nucleon. At small temperatures 
T the baryon contribution to <£> is strongly suppressed by the Boltzmann factor e~ E l T and is negligible. Below 
we assume that the proportionality m ~ 77 1 / 3 is valid in a broad temperature interval. Arguments in favor of 
such an assumption are based on the expectation that the baryon masses vanish at T = T p h(p) and on the 
dimensional grounds. Near the phase transition point E = ^/ p 2 + m 2 ~ p + c 2 rf^ /(2p). At 77 — > all baryons 
are accumulating near zero mass and a summation over all baryons gives us - instead of eq. - the following: 

*(*/, p) = $0 + \at rf + -B rf + Crf^p , (3) 

where C = Y^i c i/(^Pi)- The term pJ2iPi 1S absorbed into $0 since it is independent on the chiral condensate 77. 
The typical momenta are of the order of the temperature, pi ~ T. Thus, Eq. (O is valid in the region 77 <C T 3 . 
In the leading approximation the term C can be considered as independent on the temperature at T ~ . 

Due to the last term in Eq. © the thermodynamical potential always have a local minimum at 77 = since 
the condensate 77 is always non-negative. At small t < there also exists a local minimum at 77 > 0, which is a 
solution of the equation: 

^1 = (at + Br, 2 )r,+ ^^ = 0. (4) 

At small enough baryon density p, Eq. Q [visualized in Figure ^a)] has, in general, two roots, 771 < 770 and 
V2 > Vo, where 770 = (—at/35) 1 / 2 is the minimum of the first term in the right-hand side of Eq. The 
calculation of the second derivative d 2 &/dr] 2 shows that the second root 772 (if exists) corresponds to minimum 
of $(77) and, therefore, is a local minimum of <£>. The point 77 = 771 corresponds to a local maximum of the 
thermodynamical potential since at this point the second derivative is always non-positive. 

The thermodynamical potential $(77, p) at (fixed) non-zero baryon density p has the form plotted in Fig- 
ure IHb). At low enough temperatures (curve T\) the potential has a global minimum at 77 > and system 
resides in the chirally broken (hadron) phase. As temperature increases the minima at 77 = and at 77 = 772 > 
becomes of equal height (curve T 2 = T p h). At this point the first order phase transition to the quark phase takes 
place. At somewhat higher temperatures, T = T3 > T p h, the 77 > minimum of the potential still exist but 
$(77 = 0) < $(772). This is a mixed phase, in which the bubbles of the hadron phase may still exist. However, as 
temperature increases further, the second minimum disappears (curve T4 = T cont ). This temperature corresponds 
to a continuous transition to the pure quark phase, in which the thermodynamic potential has the form T5. 

Let us calculate the temperature of the phase transition, T p h(p), at non-zero baryon density p. The transition 
corresponds to the curve T2 in Figure db), which is defined by the equation $(772, p) = $(77 = 0,p), where 772 is 
the second root of Eq. Q as discussed above. The solution is 

/ v 3/5 / v 2/5 



and the second minimum of the thermodynamic potential is at 772 = [4a (T^ — T p h(p))/(5 B)] 
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Figure 1: (a) Graphical representation of Eq. "I" is the first term and "II" the second term (with the 
opposite sign) in the r.h.s. of the equation, (b) The thermodynamic potential vs. the chiral condensate at 
a fixed baryon density p > 0. At low enough temperatures, T = Ti, the system resides in the chirally broken 
(hadron) phase. The first order phase transition to the quark phase takes place at T p h = Tb > Ti. At somewhat 
higher temperatures, T3 > T p h the system is in a mixed state. The temperature T4 = T CO nt corresponds to a 
continuous transition to the pure quark phase, in which the thermodynamic potential has the form T5. 

At a temperature slightly higher than T p h(p) the potential is minimal at 77 = 0, but it has also an unstable 
minimum at some 77 > 0. The existence of met ast able state is also a common feature of the first order phase 
transition (e.g., the overheated liquid in case of liquid-gas system). With a further increase of the density p 
(at a given temperature) the intersection of the two curves in Figure [IJa) disappears and the two curves only 
touch one another at one point 77 = 774. At this temperature a continuous transition (crossover) takes place. The 
corresponding potential has the characteristic form denoted as T4 in Figure ^a). The temperature T4 = T cont is 
defined by the condition that the first (@J) and the second derivatives of Eq. ((3} vanish: 



and the value of the chiral condensate, where the second local minimum of the potential disappears is given by 

(o) V 2 

774 = [2a(T cont (p) — Tc )/(5B)] . At temperatures T > T cont (p) the potential has only one minimum and the 
matter is in the state with the restored chiral symmetry. Thus, in QCD with massless quarks the type of phase 
transition with the restoration of the chiral symmetry strongly depends on the value of baryonic density p. At a 
fixed temperature, T < Tj°\ the phase transition happens at a certain critical density, p p h- According to Eq. © 
the critical density has a kind of a "universal" dependence on the temperature, p p h(T) [T^ — T] 5 / 3 , the power 
of which does not depend on the parameters of the thermodynamic potential, a and B. 

The expected phase diagram is shown qualitatively in Figure|2Ja). This diagram does not contain an end-point 
which was found in lattice simulations of the QCD with a finite chemical potential ^J^j. We expect that this 
happens because in our approach a possible influence of the confinement on the order of the chiral restoration 
transition was ignored. Intuitively, it seems that at low baryon densities such influence is absent indeed: the de- 
confinement phenomenon refers to the large quark-anti-quark separations while the restoration of the of the chiral 
symmetry appears due to fluctuations of the gluonic fields in the vicinity of the quark. However, the confinement 
phenomenon dictates the value of the baryon size which can not be ignored at high baryon densities, when the 
baryons are overlapping. If the melting of the baryons happens in the hadron phase depicted in Figure Eta), 
then at high enough density the nature of the transition could be changed. This may give rise in appearance of 
the end-point observed in Ref. [El EH- The domain where the inequality \at\ ^> Cprj 2 ^, p 7^ is fulfilled, has 
specific features. In this domain the phase transition looks like a smeared second order phase transition: the 
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Figure 2: The qualitative phase diagram at finite baryon density and temperature based on the analysis (a) 
without and (b) with indication of the approximate 2-nd order transition domain. 



specific heat has (approximately) a discontinuity at the phase transition point, AC P — a 2 T c /B. The correlation 

length increases as (T - 7T ; ) at T — T^ J 0. The latter arises if we include the derivative terms in the 
effective thermodynamical potential. The phase diagram with this domain indicated may look as it is shown in 
Figure E^b). Note that the applicability of our considerations is limited to the region |T — T c (0) |/T c (0) < 1 and 
low baryon densities. 

In the real QCD the massive heavy quarks (the quarks c, 6, t) do not influence on this conclusion, since their 
concentration in the vicinity of T w Tj ^ ~ 200 MeV is small. However, the strange quarks, the mass of which 
m s ~ 150 MeV is just of order of expected Tj°\ may change the situation. This problem deserves further 
investigation. 
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